Surface Screening in the Casimir Force 
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We calculate the corrections to the Casimir force between two metals due to the spatial dispersion 
of their response functions. We employ model-independent expressions for the force in terms of the 
optical coefficients. We express the non-local corrections to the Fresnel coefficients employing the 
surface d± parameter, which accounts for the distribution of the surface screening charge. Within 
a self-consistent jellium calculation, spatial dispersion increases the Casimir force significatively for 
small separations. The nonlocal correction has the opposite sign than previously predicted employing 
hydrodynamic models and assuming abruptly terminated surfaces. 
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The Casimir force between two ideal mirrors a small distance apart is a manifestation of the quantum fluctuations of 
the electromagnetic field confined to the cavity within pj . The study of vacuum forces among real materials began with 
Lifshitz 0] in 1956, who obtained a formula for the force between semi-infinite homogeneous isotropic local systems 
characterized by a frequency lu dependent dielectric function e(w). Developments in micro- and nano-technologies 
have stimulated a renewed attention on the Casimir force. Deviations from the force between ideal mirrors are now 
measurable experimentally. The pioneering measurements of Lamoreux .33 bad a 5% precision and 1% precision has 
now become common 0,|a,IH0- Distances down to « 60nm have been explored Due to alignment difficulties, 
most experiments have involved spherical surfaces, while most theories have concentrated on the force between flat 
parallel surfaces for simplicity, so comparisons have resorted to the P ro Xl rmty theorem Q which relates approximately 
the force between curved surfaces to that between flat parallel surfaces. Nevertheless, measurements of the force 
between flat surfaces have already been performed [ToT ] . We expect that in the near future even better measurements 
at smaller distances will be produced. Therefore, manifold corrections have to be carefully evaluated in theoretical 
calculations, such as the finite conductivity of the metals [III Ha . Ha . . surface roughness |15l[l6j . finite temperatures 
Ell [13 grain structure 0, etc. In this paper we calculate non-local effects in the Casimir force. We employ density 
functional jellium calculations [T^ . l2pj of the surface response functions which are consistent with realistic electronic 
density profiles that decay smoothly to zero beyond the nominal surface of the metal, instead of being artificially 
truncated abruptly |2lll22 |. At very small distances the non-local corrections can be remarkably large. Within simple 
hydrodynamic models with abrupt surfaces |2ll |22| nonlocal effects yield a diminution of the force apparently due to 
the possibility of exciting bulk plasmons at the surface of spatially dispersive metals. However, nonlocal corrections 
within self-consistent jellium theories have the opposite sign and yield an increase in the Casimir force, as most of the 
surface screeening charge accumulates outside of the body of the conductor, as characterized with d± theory list l20l . 

It has been shown [2lL l23l ] that Lifshitz formula 0] , when written in terms of the exact surface impedances )2 li l25j | 
of the system, or equivalently, in terms of its optical reflection amplitudes |26| . 
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is applicable to any system with translational invariance along the surfaces and isotropy around their normal |21l l'2.'"{l| . 
Here, F(L)/A is the force per unit area between two slabs with flat parallel faces a distance L apart, c the speed 
of light, q = lo/c the wavenumber, Q the projection of the wavevector parallel to the surface and k = \J ' q 2 — Q 2 its 
component normal to the surface, £ Q = (rl l r 2 Y e 2lkL )~ 1 , rJJ the reflection amplitude of slab a — 1,2 corresponding to 
the polarization a = s,p, and / = coth(/3?kj/2)/2 the photon occupation number at temperature T = l/fcs/3. The 
integral over k in Eq. goes from iQ towards along the imaginary axis and then towards oo just above the positive 
real axis, although it may be easily manipulated into a more convenient path along the imaginary axis. 

Nonlocal effects l25| in the Casimir force have been calculated approximately 1271 employing expressions for the 
surface impedance [l3| which are only valid for good conductors at low frequencies |28j |. and they have been calculated 
from the density response function of semiinfinite metals within a semiclassical inifinite barrier model |22^ . However, 
they may be fully incorporated simply by plugging the appropriate reflection amplitudes into Eq. . The polarization 
induced at a given position r within a non-local or spatially dispersive medium depends on the field that is applied 
at nearby positions f* [2{| , instead of depending exclusively on the field at f. Within the bulk of a homogeneous 
system, non-locality manifests itself as a dependence of the response functions on the wavevector besides their usual 
dependence on the frequency. Thus, spatial dispersion is expected to become important as the lengthscales associated 
to the field become small. A few consequences of spatial dispersion at metal surfaces, such as the anomalous skin effect 
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[30| , the excitation at surfaces of propagating plasmons [21J , and the finite screening distance |23 of low frequency 
fields, have been studied in relation to the Casimir force, which they modify substantially at very small distances. 

The most simple non-local model for the response of a metal is the hydrodynamic model m which a negative 

semiclassical electronic fluid is confined within a uniform positively charged region. The electrons interact with the 
macroscopic electric field and with the hydrodynamic forces derived from the pressure of the fluid. The longitudinal 
contribution to the bulk dielectric response in this model is 
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where the compressibility, related to (3 2 — 3vjp/5, is the source of the spatial dispersion. Here, vf is Fermi's velocity, ui v 
the plasma frequency, r a phenomenological lifetime, and k\ is the component of the wavevector of the longitudinal wave 
normal to the surface, given by the longitudinal dispersion relation e L = 0. In this model, the reflection amplitude 
r s for s polarized light agrees with that, r°, of a local metal described by the Drude dielectric funcion e T (uj) = 
1 — lu 2 /(lo 2 +iw/r). To calculate the reflection amplitude r p for p polarized waves, the excitation of longitudinal 
waves at the surface has to be accounted for. This might be done by postulating the continuity of the component Ej_ 
of the electric field normal to the surface as a physically reasonable additional boundary condition (meaning there 
is no singularity in the induced charge density), beyond the usual conditions derived from electromagnetic theory, 
and solving them for the amplitudes of the longitudinal transmitted field as well as the transmitted and reflected 
transverse waves, yielding 
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where k m = \J e T q 2 — Q 2 is the normal component of the wavevector of the transverse wave transmitted into the 
metal 33] . The reflection amplitude r° for a local metal may be obtained from Eq. © by taking the limit f3 — ► 0, 
ki — > oo. 

In order to use more realistic models of metallic surfaces, we require expressions for the reflection amplitudes 
that may also be employed beyond the hydrodynamic model. We remark that non-local effects are expected to be 
important mostly within a thin region close to the surface of the metal, where the electric field varies rapidly from 
its value in vacuum to its bulk value. As the width of the selvedge is typically much smaller than the wavelength A, 
r p may be calculated perturbativelv l34l . Applying boundary conditions across the selvedge and employing the long 
wavelength approximation we obtain \Hk |20| 
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where d±/X is the small perturbative parameter and 
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is the position of the centroid of the induced charge density 5p(z) measured towards vacuum from the nominal surface 
of the metal defined by the edge of the positive background. In a local model with an abrupt interface, there is a 
singular charge induced right at the surface, so d± = 0. In non-local models with an abrupt termination the screening 
charge occupies a finite region within the metal, so that d± < 0, as may be confirmed within the hydrodynamic model, 
for which 
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In this case, — d± may be regarded as the finite frequency generalization of the Thomas-Fermi screening distance. 
However, within more sophisticated models such as the self consistent jellium model \\s\ . l2(| with a smoothly termi- 
nated charge density profile, the screening charge may be localized on the outside of the surface, so d± may become 
positive. Furthermore, surface dissipation yields a complex valued d±. 

Linearizing the Casimir force Q with respect to d± , we obtain an expression for the non-local correction 
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FIG. 1: Nonlocal correction to the Casimr force between Au flat surfaces (a-d) and between a flat and a curved surface (e), 
calculated within the hydrodynamic model using the exact (a) reflection amplitudes, retarded d-parameter theory using the 
dynamical values of d±(u>) (b, e) and its static limit d±(0) (c), and the non-retarded d-parameter theory (d) as a function of the 
separation L between plates. The left panel shows the normalized correction —SF/F and the right panel shows the correction 
(—6F/F)(Lui p /c) scaled by the separation L. 



where £° is the local limit of £ p . This expression may be simplified at small distances for which retardation may be 
neglected, yielding 

sf(l) 2h C r, nn3 r, „ 

where we may further approximate = [(e + l)/(e — l)] 2 e 2< 2 L . 

Fig. QJ, calculated within the hydrodynamic model using the density parameter r s — 1.59A and the lifetime 
t = 400/wp corresponding to Auj^a], shows that the non-local correction decreases the Casimir force within the 
hydrodynamic model. At L — 60nm rs 3c/cj p , the closest distance for which measurements are available, the non- 
local correction amounts to 8F/F rs 0.5%, which is smaller |33( than the experimental error estimate of 1.75% |8(. 
However, 5F/F increases for shorter distances approximately as and may become of the order of 100% at a few 
A. Although non-local corrections seem too small to be noticeable in current experiments, they will become much 
more important at closer range. 

We use now the hydrodynamic model as a test ground to appraise the results of the perturbative theory In order 
to enhance the differences between the different approximations discussed above, in the right panel of Fig. Q]we show 
the corresponding nonlocal corrections scaled by the separation L. For distances larger than L rs 2nm ri 0.1c/w p 
there is a very good agreement between the perturbative and the exact hydrodynamic calculations; for that distance 
perturbative calculation produces a Ri 12% error on the rs 15% nonlocal correction. We also show in Fig. [I] that the 
corrections to the force between a plane and a sphere |8| are similar to those for flat surfaces. 

The non retarded calculation also yields similarly good results up to distances L rs 20nm rs c/u> p , and although 
it overestimates the corrections by up to rs 20% at larger distances L rs 10c/cj p , they are themselves negligible 
SF/F < 0.1% at those distances. Finally, we calculated the nonlocal corrections employing the perturbative theory 
but substituting the static value d±(0) instead of the dynamic value d±(u>). In this case, the nonlocal correction is 
further overestimated but by not more than 20% if L > {).lc/uo p . 

As the vacuum fluctuations of the electromagnetic field within a cavity are correlated to fluctuations in the charges 
and currents induced at the surfaces, the p-polarization contribution to the Casimir force may be interpreted as due 
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FIG. 2: Normalized nonlocal corrections to the Casimir force calculated for the selfconsistent semiinfinite jellium model using 
the LDA for different densities corresponding to r a /clb = 2, 3, 4, and 5. 



to the force between fluctuating charges in one surface and their images charges at the opposite surface. Within 
nonlocal optics, the charges induced at a surface are not at the nominal surface, but are actually spread over a finite 
region whose centroid lies a distance d± from the nominal surface. At low frequencies, the position of the centroid 
corresponds to the effective image plane |19|. Thus, we may expect that the force between nonlocal media is similar 
to the force between local media but separated by an effective distance V = L + SL. As F oc L~ 3 at small distances, 
SF/F rs —35L/L, which allows to identify 5L w — 0.006c/o-> p sw — 2d±(0) from Fig. ^ Thus, nonlocal effects amount 
to displacing each surface a distance ~ d± towards each other. Since in the hydrodynamic model d± < 0, the effective 
separation is larger than the nominal separation and nonlocality decreases the force. 

Although very simple and amenable to analytical solutions, the hydrodynamic model for an abruptly terminated 
electron gas is not considered a realistic model of metal surfaces. It doesn't account for the fact that electrons spill 
across the nominal surface, giving rise to a static surface dipole which is the source of the confinement potential, nor the 
density oscillations due to the quantum interference between the incoming and outgoing wavefunctions of electrons 
reflected at the surface. Landau damping due to excitation of electron-hole pairs and many body effects such as 
exchange and correlation are also absent. Nevertheless, much more sophisticated models of metal surfaces have been 
developed [l9ll20) . and have been employed to calculate d± for manifold metals. In the jellium model, the Schrodinger- 
like Kohn-Sham equations [3|J of density functional theory are solved within the local density approximation (LDA) 
to obtain the selfconsistent wavefunctions and ground state density of electrons which interact with an abruptly 
terminated semiinfinite homogeneous positive background. Their response to a perturbing electric field may then be 
obtained using time dependent LDA, and d± may be calculated from the ensuing density-density response function 
[U using Eq. ©. 

In Fig. J2)l we show the nonlocal corrections to the Casimir force for metals of various density parameters r s 133 
calculated within the self consistent jellium model obtained by plugging the corresponding known values of d±(0) [2(j 
into Eq. J7J. In contrast to the previous hydrodynamic results, in this model the nonlocal effects increase the Casimir 
force since most of the screening takes place in the low density region on the outside of the nominal surface, i.e., 
d± > and the effective distance L' is smaller than the nominal distance L. We remark that in the hydrodynamic 
model the ground state electronic density is constant up to the surface, where it is abruptly truncated, while in the 
jellium model it has a tail that extends beyond the surface. The stiffness of the tail of the electronic fluid is smaller 
than within the metal, so that it is more easily polarized than the bulk. Therefore, the centroid of the induced charge 
within the jellium model lies outside of the metal. In the hydrodynamic model there is no charge at all out of the 
metal, so the centroid of the induced charge is necessarily within the metal and has the opposite sign. According to 
Fig. [3 the nonlocal corrections are larger for high density metals (small r s ) than for low density metals (high r s ). 

In summary, we have calculated the nonlocal corrections to the Casimir force between metals using a generalization 
of Lifshitz formula in terms of reflection amplitudes which has been proved to be applicable to arbitrary materials 
that are homogeneous and isotropic along the surface. A perturbative theory in which the width of the selvedge 
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divided by the relevant wavelengths plays the role of small parameter is a useful approximation for distances that are 
not too small, and yields an interpretation of the nonlocal effect in terms of an effective separation between the metal 
surfaces. The exactly solvable hydrodynamic model predicts a diminution in the force as the effective separation 
increases due to the finite screening length of spatially dispersive metals. Contrariwise, the jellium model yields an 
increase of the Casimir force due to the spilling of electrons out of the metals and into vacuum, and to the fact that 
most of the screening takes place in the resulting low density region outside of the nominal surfaces, thus reducing 
the effective separation. \5F/F\ grows roughly as L _1 and thus becomes very important at distances smaller than 
~ fOOnm. As the surface screening may be manipulated through the adsorption of overlayers 20], our results suggest 
that the Casimir force at small distances may be tailored. Other surface corrections, such as the surface local field 
effect [37j may easily be incorporated into Casimir force calculations simply by replacing the Fresnel coefficients r p 
and r s in Eq. by the appropriate nonlocal, surface-corrected expressions [34J instead of Eqs. © or Q. We hope 
our results stimulate experiments exploring such small distances. 

We acknowledge useful discussions with Carlos Villarreal, Raul Esquivel, Gregorio H.-Cocoletzi, Catalina Lopez- 
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